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DISCUSSION OF EARTHQUAKE STRESSES IN SHEAR BUILDINGS 
PROCEEDINGS-SEPARATE 177 


Joun A. Biume,'* M. ASCE.—There has been considerable technical litera- 

16 Cons. Civ. and Structural Engr., San Francisco, Calif 
ture concerning the analysis of buildings subjected to earthquakes. A great 
deal of this revived interest in building dynamics has undoubtedly been induced 
by the possible application of earthquake-resistant design methods to blast de- 
sign. Such renewed interest is heartily welcomed. A word of caution is re- 
quired, however, regarding the general acceptance of mathematical or numerical 
methods based on idealized structures as constituting complete earthquake- 
resistant design. Caution is especially necessary when the buildings are con- 
veniently subjected, for the purpose of analysis, to sinusoidal motion of specific 
characteristics, or even to a few cycles of sinusoidal motion of different char- 
acteristics. 

The author’s contribution can be considered another method of investigat- 
ing building dynamics, but neither this paper, that by C.S. Whitney, M. ASCE, 
B. G. Anderson, and the author,'® nor any other idealized mathematical treat- 


"4 Comprehensive Numerical Method for the Analysis of Earthquake Resisting Structured,’ by 
C. 8. Whitney, B. G. Anderson, and M. G. Salvadori, Journal, A.C_1., Vol. 23, No. 1, September, 1951. 


ment should be thought to constitute a panacea. In this regard, reference is 
made to the statement in Mr. Salvadori’s conclusions: 


“The availability of such procedures should encourage the designer to 
abandon the simplified and unrealistic methods of analvsis suggested by 


codes (such as the 0.1g method) with obvious gains in economy and safety.” 


The mere assumption of static lateral forces, whether based on 0.1g or any 
other percentage of g, is an unrealistic approach to a dynamic problem. Build- 
ing codes, including the San Francisco (Calif.) code and the 1952 Uniform 
Building Code,"’ vary the seismic factor with the height or number of stories. 


Uniform Building Code,’ Pacific Coast Building Officials Conference, Los Angeles, Calif., Vol. I, 
1952, p. 280. 


A suggested code'* approaches the problem from a more realistic viewpoint— 


18TLateral Forces of Earthquake and Wind,” Transactions, ASE, Vol. 117, 1952, Appendix I.— 
Lateral Force Code, p. 746. 


permitting decreases in lateral design forces with flexibility and considering 
actual damage and earthquake records. It is untrue that the use of assumed 
static lateral forces in design is entirely or evea largely caused by the lack of 
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numerical procedures for dvnamic analvsis. There are actually many methods 
and techniques—mathematical, experimental, and empirical. Two math- 
ematical methods which can be used to obtain building periods in much less 
time than the method suggested by Mr. Salvadori will be presented. 

Some of the factors involved in the design of buildings to resist earthquakes 
greatly complicate and make questionable the implied accuracy and security 
of any rigorous analysis as a sole criterion. These factors are: 


1. The motion of the earth in the general area of the greatest energy release 
is not a simple sinusoidal movement or even a superimposition of only a few 
such waves of variable period, amplitude, and acceleration. The motion is 
violent, chaotic, and in three directions. It may last from 30 sec to 40 sec or 
more. The motion is exceedingly complex and variable. 

2. The epicenter is not known in advance, although active faults usually are. 

3. There are many active faults that may affect structures in some locali- 
ties. The assumption, therefore, that distant shocks may transmit smoother 
waves of lesser intensity (but more subject to analysis) can be dangerous as 
the only basis for design. 

4. The natural periods of the building, which are greatly affected by the 
walls, fireproofing, ground conditions, and adjacent structures, are subject to 
change during an earthquake. With increased amplitude the periods shorten 
as more partitions, walls, and other items come into structural action, and the 
periods lengthen as cracks and failure occur. 

5. The damping is not only conveniently viscous but is frictional and varies 
not only with materials and types of construction but with amplitude and 
damage. 

6. The concept of using dominant ground periods or important Fourier 
components of actual records as disturbing forces or motions is subject to more 
exploration before being a reliable basis for design. Moreover, such character- 
isties may be not only local, but local for a specific origin of energy release, 
which point may not be the same in future shocks. 


There are several other factors complicating rigorous aseismic design. This 
approach should be joined with (a) the study of actual damage (and lack of 
damage) in exposed structures; (b) field research; (c) dynamic laboratory re- 
search (shaking table and similar equipment); (d) mathematical research; (e) 
obtaining and analyzing reliable earthquake records, particularly for the re- 
sponse of various systems of one and more masses to such motion, with various 
damping coefficients, and spectrum studies; (f) practical design considerations ; 
(x) practical considerations of function, economy, and risk; (h) complete 
drawings showing details; {i) adequate construction supervision; (j) code re- 
quirements; and (k) committee work and discussion. 

The writer has been engaged in all of these phases of the problem and wishes 
to defend (at- the present state of knowledge) modern (1953) methods using 
assumed lateral forces for codes. It is suggested that such codes be constantly 
improved as experience and knowledge indicate. The static lateral forces can 
be, and to some extent have been, based on dynamics and mathematics as well 
as on many other considerations. Special structures are often designed dy- 
namically. Earthquake-resistant design is essentially an art which uses either 
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assumed static forces or dynamic forces (or both) as a means to an end. No 
code, mathematical method, or assumed force can remove the basic requirement 
for sound, practical judgment and experience. Above all, continued research 
is needed. 

Research was performed at Stanford University at Stanford, Calif., in 1933 
with Lydik 8. Jacobsen as faculty sponsor."* This work disclosed many inter- 

'* The Reconciliation of the Computed and Observed Perinds of Vibration of a Fifteen-Story Build- 
ing,” by John August Blume and Harrv Leonard Hesselinever, thesis presented to Stanford University, at 
Stanford, Calif., in 1934, in partial fulfilment of the requirements for the degree of Engineer. 
esting factors concerning the dynamic characteristics of a building and the 
methods of computing these characteristics. The results of this research were 
not widely distributed because they were considered interesting but academic 
in view of the many other factors in the design problem. Many of the writer’s 
comments on the author’s paper are based on that research. 

Shear versus Flerure.—The natural vibration frequencies of a hypothetical 
five-story structure have heen developed by Mr. Salvadori from generalized 
equations based on certain assumptions. There is little error involved in the 
assumption that the mass of each story is concentrated at the corresponding 
floor level when four or more stories are involved. This assumption was the 
construction basis of dynamic building models at Stanford University by Mr. 
Jacobsen. However, in the introduction the author states that a structure with 
an height-to-width ratio of 5 or less ean he termed a “‘shear building,” and the 
influence of direct stresses (flexure of the building as a whole) can be neglected. 
The listed reference? does not reconcile the difference between wind and earth- 

2"*Wind Bracing.’ by H. V. Spurr, McGraw-Hill Book Co., New York, N. Y. 1930. 
quake design. Normal wind design for framed buildings considers frame re- 
sistance to all wind forces, whereas earthquake design usually considers not 
only the frame, but fireproofing, walls, and other structural elements—in pro- 
portion to their relative rigidities. Of course, there are some exceptions in 
both cases. If the walls and fireproofing are neglected, the shear distortion 
between level floors is much greater than the deflection caused by direct stress 
in the columns, and the 5-to-1 ratio applies. In earthquake design, however, 
the walls, floors, fireproofing, permanent partitions, and shear walls provide 
not only mass but rigidity and strength. Many practicing engineers are con- 
vinced of the value of a good frame as a secondary defense,'* but they do not 
neglect the walls or other structural elements which can provide both strength 
and rigidity and which may constitute a hazard to personnel. 

An actual fifteen-story building was analvzed in detail at Stanford Uni- 
versity."* The height-to-width ratio in one direction was 3.27. The funda- 
mental period computed for this same direction with shear alone (neglecting 
flexure and ground movement) was 0.605 sec as compared to 1.026 sec if shear 
and flexure were both considered. The error of the former is 70%. When 
ground flexibility is introduced the period becomes 1.32 sec. In the other 
direction, the height-to-width ratio was 2.86 and the fundamental period with 
shear alone was 0.764 sec as compared to 1.04 sec with shear and flexure; the 
error was therefore 36%. With ground movement in this direction, T; was 
1.23 sec. An error in T, for the same building caused by neglecting flexure 
would be from 15% to 20%. The walls of this building contained a consider- 
able number of windows—about 47% of the wall length in the first direction 
and 60% in the other direction. 
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ratios, flexure as a whole should be considered in the determination of the first, 


For buildings with an height-to-width ratio of 5, and with other slenderness 


second, and possibly the third modes. 

The First Mode.—The first mode of an hypothetical building is evaluated in 
Table 4 by the Stodola-Vianello method. Other methods give results with less 
expenditure of work and time and with greater ease of introducing factors such 
as flexure, ground rotation, and ground translation. One simple solution for 
the first mode for any number of different masses and for different spring con- 
stants can be made directly without trial and with only elementary mathe- 
matics; the Lord Rayleigh potential energy method is applied to d’Alembert’s 
principle and the basic equations of motion. Although mathematically approx- 
imate, the solution can be considered sufficiently accurate for engineering pur- 
poses, which involve many factors for which the value of the second significant 
figure is doubtful. 

This method will be demonstrated by using the same geometric and elastic 
properties as set forth in Table 1. It is assumed that the dynamic deflection 
curve for this mode is the same as the static deflection curve of the system 
loaded laterally by its own weight. Table 9 illustrates the entire process and 


TABLE 9.—Direct CompuTATION OF APPROXIMATE 
First-Mope FREQUENCY 


‘ |W =Miol x da, = at W at, 


(10) 
323,604 


(7) (8) 


5 11.65 8.47 71.88 38, 

4 15.53 6,001 10,503 7.98 1.316 7.820 61.15 46,928 366,961 
3 15.53 6,001 16,504 9.12 1.810 6.504 42.30 39,031 253,842 
2 15.53 6,001 22,505 10.26 2.193 4.694 22.03 28,169 132,202 
1 15.53 1 28,506 11.40 2.501 2.501 6.26 15,009 37,506 


167,305 1,114,115 


has an extra column, M,, which would ordinarily be eliminated since W would 
be entered directly in pounds or kips. The equation, 


can easily be proved valid by equating the maximum potential energy in the 
fully deflected position to the maximum kinetic energy when vibrating through 
the static equilibrium position. It is simple to add the flexural deflection, Ay, 
and ground rotation values, 4,, to the shear deflecticns before computing W A 
and W 

From Table 9 and Eq. 71 


167,305 


w*, = 386.4 X 1,114,115 


from which 
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The value of w, computed in one step is 0.79% in error from the value (7.557) 
determined from Eq. 5. This error from the classical solution is much less than 
the knowledge of the actual modulus of elasticity of concrete or steel. 

It should be noted that w, (‘frequency” in the author’s nomenclature) is 


actually the angular velocity in radians per second. The frequency is =. 


The first period, 7;, is 2m which is 0.824 sec according to results of Table 9. 


Fig. 3 is a qualitative plot of the first-mode curve, wherein the values 
obtained from Table 9 are reduced in proportion to the maxima obtained by the 
two methods. 

A more accurate method of computing w, and all natural frequencies involves 
a trial type of solution, but less labor than Mr. Salvadori’s procedure. This 
procedure’ will be demonstrated subsequently. It was adapted to building 

Die Berechnung der Drehschwingungen,"’ by H. Holzer, J. Springer Verlag, Berlin, 1921. 
dynamics (apparently for the first time) from a device first given by H. Holzer™ 
and applied to the torsional vibration problem of a shaft by 8. Timoshenko."'. 


Vibration Problems in Engineering,’’ by 8. Timoshenko, D. Van Nostrand Co., Inc., New York, 
N. Y., 2d Ed. 1937, p. 263. 


General Method for Shear Vibration.—For convenience and flexibility in 
introducing variable ground characteristics in the research at Stanford Uni- 
versity,"® the top story (in lieu of the lowest story) was called | in the i series 
and the nth story the lowest. In order to avoid confusion with the author's 
notation, the notation z is used instead of i (see Fig. 4). 

From d’Alembert’s principle, 


M,2,+ K, (2. — 2141) — (22-1 — 2,) = 0 


and by the addition of the n equations of motion. 
M,z, + Kz (tn) = 0 


Assuming simple harmonic motion (for the natural system, not the forcing 
distrubance), 
X,sinwl....., 
and 
Z,=— X,w'sinwlt 


Substituting and letting sin w ¢ equal unity, for maximum motion, 


= 0 


Also, for the zth mass it can be shown that 


1 e-1 


X, = X..1- 


| 
n 
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The process can best be undertaken in tabular form as in Table 10, which 
uses the same geometric and elastic values established by the author. 


TABLE 10.—Examp.Le or Frequency CALCULATION WITH w 
AssuMED Equa. To 20.61 


x 


4948.6 

1824.0 

—3774.0 

—5941.6 

—0.6138 —4049.1 


The first step in this process is to assume an w-value and then to assume 
that X is equal to unity. Since frequencies of such a system are independent 
of amplitude, the X-values can later be proportionally adjusted, if desired. 
At a solution for an harmonic period, Eq. 77 must be satisfied. 


6992.1 + 11,400 (—0.6138) 5 = 0 


In this case (in Table 10) the solution is for the second mode since the criterion 
is essentially met and the X-values indicate two nodes, or second-mode move- 
ment. The procedure is to select various values of w and to plot the results. 
With a few trials, which take little time with a calculating machine, solutions 
for frequency as well as mode shapes are readily obtained. 

The curve in Fig. 5 is the curve of Eq. 77 for the author’s building. It is 
ordinarily unnecessary, of course, to plot the curve except near the desired 


frequency. 

It should be noted in Fig. 5 that each crossing of the vertical zero axis 
represents a natural-mode frequency. All those shown agree with those com- 
puted by Mr. Salvadori after a more laborious process. In Fig. 3 the shape of 
the fundamental mode coincides exactly with the curve obtained from Eq. 6b. 
The X-values of Table 10, when proportioned to the author’s top-story value, 
agree with his second-mode values of Table 3(a). 

The mathematical research undertaken at Stanford]University did not limit 
the computation of the first or second modes to shear only, but included flexure 
of the building as a whole, ground rotation, ground translation, the effectiveness 
of brick filler walls in conjunction with the steel frame, and other considerations. 
Actually the computation methods were means to the end of learning how{multi- 
story buildings vibrate. It has been most interesting to substitute Mr. Salva- 
dori’s hypothetical building values in the equations and in the methods de- 
veloped in 1933, and to obtain the same results as did the author. 


3 
ed (1) (2) (3) (4) (5) (6) (7) (8) 
f 5 11.65 4948.6 4948.6 0.1462 0.7235 
4 15.53 6506.7 6772.6 0.1253 0.8486 
s 3 15.53 6596.7 2098 6. 0.1006 0.3286 
2 15.53 6596.7 —2943.0 0.0975 —0.2869 
1 15.53 6596.7 —6992.1 0.0877 
~ 
Fe. 
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Fia. 5 


E. F. Masur.”—It is gratifying to note that Mr. Salvadori’s paper marks 

™ Associate Prof. of Civ. Eng., Illinois Inst. of Technology, Chicago, II. 

a return to the “classical” procedure of determining the dynamic response of an 
elustic system to a set of applied forces by expanding that response in terms of 
the natural modes of vibration of the system. This approach has the further 
advantage (in addition to the ease of computation mentioned by the author) of 
permitting the study of various loading conditions with a minimum of addi- 
tional effort. 

It is felt that the concept of a shear building needs further elucidation. 
Other than the two conditions stated in the Introduction, the author also as- 
sumes that the floors are infinitely rigid. This simplification leads to the almost 
diagonal form of the frequency determinant in Eq. 5. There is some doubt, 
however, whether the error thus introduced is negligible in the case of actual 
structural-design practice. In fact, the existence of joint rotations lessens the 
resistance of a structure against sidesway to such an extent that it is reasonable 
to expect appreciably smaller natural frequencies of vibration than those based 
on the idealization of complete floor rigidity. This weakening effect is well 
recognized in conventional wind-bracing analysis, in which the stiffness of the 
floor girders is as much a factor as that of the columns. 

To illustrate this point numerically the simple one-bay, three-story frame 
shown in Fig. 6 should be considered. The mass, modulus of elasticity, and 
moment of inertia are, respectively, M, E, and J for the columns, and 8 M, E, 
and 4 / for the girders. When the distributed spandrel masses are combined 
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with the respective floor masses, the frequency equation, on the assumption 
of rigid girders, takes the form: 


(48 D — 10 w*) — 24D 0 
— 24D 48 D — 10a? — 24D)| =0.... (79a) 
0 24D 


in which 


The roots of Eq. 79a—that is, the natural frequencies of the idealized structure, 
are 


and 


These results differ from the exact values (that is, “exact’’ within the other 
two assumptions regarding shear buildings) of 


and 


It can be seen that the nonrigidity of the floors affects the magnitude of the 
natural frequencies of vibration to a significant extent. Furthermore, it can be 
shown that the effect on the modes is as great, or greater, than the effect on the 
natural frequencies of vibration. 

It is possible, however, to make some allowance for the joint rotations with- 
out sacrificing the simplicity of the form of the author’s equations of motion. 
This can be done by assuming that all columns and girders display a point of 
contraflexure at their midpoints. In that case, the conventional analysis 
shows that the spring constant A,, as used by Mr. Salvadori, must be multi- 
plied by a factor, 


in which >} ¢ extends over all the girders and >> y extends over all the members 
(that is, girders and columns) meeting at a joint. In Eq. 82, S is used to 
denote the stiffness of a member. 

If this procedure is applied to the example under consideration, the fre- 
quency equation is modified to 


(30 D — 10) - 12D 0 | 
-— 12D (26 D — 10 w*) -14D =0 
0 - 14D 14D 
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The roots of Eq. 83 are 


and 
= 4.37 D 


It is apparent that the proposed method vields a good approximation of the 
fundamehtal frequency (and, to a lesser extent, of the fundamental mode) of 
vibration. Conversely, the accuracy of the method at the high frequencies is 
less than the accuracy obtained by the author. This lessening of the accuracy 
appears reasonable if it is considered that the higher modes involve nodes which 
reduce the magnitude of the joint rotations. 


= 
~ 
bes Fia. 6 
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Asst. Engr., Ammann & Whitney, New York, 
™ Structural Designer, Ammann & Whitney, New York, N. Y. 


Leon 8S. Levy.**—A valuable contribution to the solution of dynamic problems 
% Professional Assistant (Mathematics), Ammann & Whitney, New York, N. Y. 
involving shear buildings has been made by the author. The writers have uti- 
lized the material presented by Mr. Salvadori for the dynamic analysis of 
structures subjected to impulsive loads. In the course of this work it was noted 
that the rotational mode, as defined by Eq. 36, is valid only for buildings with 
large height-to-width ratios where the distribution of mass across the width can 
be neglected. However, under these conditions, the number of actual modes and 
frequencies in a building subject to shear and rocking will not exceed the num- 
ber of stories. As described by Mr. Bleich,?* the 
%* Frequency Analysis of Beam and Girder Floors,” by Hans I. Bleich, Transactions, ASCE, Vol. 115, 


1950, p. 1035 


**e** total number [of frequencies] can include a trivial root w = =, 
however [and a mode equal to zero), which happens in the special case in 
which the complete system contains only masses which already form part 
of system B.” 


For structures in which the distribution of mass across the width cannot be 
neglected, the actual rotational modes can be derived as follows (in the same 
manner as the shear modes): 

The derivation of the modes from Eq. 33, 


16+ B6=0 


will be valid only when the proper dimensionality has been introduced.”’ 


” “Mathematical Methods in Engineering,’ by Tn. von Karman and M. A. Biot, McGraw-Hill Book 
Co., Inc., New York, N. Y., 1940, p. 177. 


Thus, setting z, equal to L, 6 and dividing Eq. 85 by L,, there results 


+ = = 1,2,---, 


Eq. 86 on application of the normalizing condition for the rotational mode, 
leads to 


Lex) Loew 


from which 


Remembering that the physical condition necessitates that r(i) = ¢ Ly, 


L, 


‘ L; 
sm. (1.42 
= ( + 
Thus, the rotational modes include a factor of the width of the building, and, 


as anticipated, the use of this factor reduces the rotational mode and the rota- 
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tion. In the illustrative example, where the building has a height-to-width 
ratio of 1, an error of approximately 15% is introduced by the use of Eq. 36. 

It should be noted that modification of the normalizing condition (following 
Eq. 416) and of the orthogonality condition also becomes necessary. 

The writers question the validity of Mr. Salvadori’s statement that 


“Inasmuch as the largest stresses due to earthquakes occur, in general, a 
short time after the beginning of the earthquake, and since the logarithmic 
decrement 7, is nearly 1, the influence of damping on elastic earthquake 
stresses is often negligible.” 


For any given earthquake the maximum structural distortions occur as the 
result of either the severe discontinuous shocks relatively close to the epicenter 
or the sustained cscillations of lower acceleration and longer period which occur 
farther away from the epicenter. 

In the case of the sustained oscillations, it is apparent that damping, even 
of small magnitude, must play an important part in reducing the effects of 
resonance for structures whose natural period approximates the period of the 
ground motions. However, even in the former case, damping cannot be ne- 
glected. It was found” that, in the case of the earthquake of May 18, 1940, at 


% Discussion by R. R. Martel, G. W Housner, and J. L.. Alford, of ‘Lateral Forces of Earthquake and 
Wind,” Transactions, ASCE, Vol. 117, 1952, p. 769, Fig. 17. 


El Centro, Calif., the response of an idealized structure could be reduced by 
more than 50% by the adaition of 2% of critical damping. Damping in ordi- 
nary buildings has been found to range from 7.5% to 14% of critical damping 
for concrete structures and 1.5% to 5% for steel structures. 

Mario G. Satvapori,” A.M. ASCE.—Mr. Blume justifiably warns the 


* Prof. of Civ. Eng., Columbia Univ., New York, N. Y., and Cons. Engr., Ammann & Whitney, 
New York, N. Y. 


engineer about the complexities of earthquake-resistant design, and the writer 
wholeheartedly agrees with him. Mr. Blume has performed a valuable service 
to the profession by listing (in factors 1 to 6 and (a) to (k)) some of the essential 
points governing the dynamics of earthquake motion in buildings. However, 
one must realize that the research engineer must confine his studies to the 
influence of some of these factors lest he be hampered by the complexity of 
the phenomenon. It is not a matter of controversy that the improvements 
in procedures stipulated in some building codes depend essentially on the results 
achieved by mathematical and experimental methods. The faster these results 
are incorporated into the codes, provided this is done intelligently, the better 
practical earthquake design will become. 

The writer’s procedures combine the classical approach with the numerical 
approach to obtain an increase in efficiency. The methods of Lord Rayleigh 
and H. Holzer are well known, but it is often a matter of familiarity that decides 
the designer’s preference. It was gratifying to learn that Mr. Blume was 
able to reproduce, by other methods, the writer’s numerical results. 

The influence of bending on the frequencies and modes of a building, which 
acts primarily as a shear building, can obviously be substantial, as shown by 
Mr. Blume. This influence can be accounted for in the procedures suggested 
by the writer. 

The writer is aware of the existence of other numerical and mathematical 
procedures which are well adapted to the solution of the earthquake problem 
and did not mean to imply that empirical methods should be disearded entirely. 
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His intention was to compare the efficiency of three methods of solving a 
simplified earthquake problem, in the hope of adding some knowledge to a 
complex field of investigation. 

Mr. Masur’s suggestion to account approximately for the influence of joint 
rotation seems to be applicable only in connection with the first mode. Since 
this mode can easily be obtained by Rayleigh’s method, the writer does not 
feel that, in general, its use would be justified. 

Messrs. Cohen, Laing, and Levy take issue with the writer’s statement 
concerning the influence of damping on earthquake motion. The influence 
of damping on the over-all response of a building depends essentially on (1) 
whether there is resonance between the Fourier components of the earthquake 
and the natural modes of the building, and (2) the amount of energy dissipated 
through friction. Although it is quite possible to encounter cases in which 
these two factors are of fundamental importance, it is not possible to state 
beforehand that this should always be the case. It would be of interest to 
study the influence of friction on earthquake motion from a mathematical 
viewpoint. 

Even though the influence of the width of the building is not among the 
most important factors in dynamic design, it can be easily taken into account 
by substituting for Eqs. 6a, 7, 36, 37, and 38, respectively, the following, in 
which @, is the rotational component of the kth mode: 


(i) Mi + > 


N 


9, = Cr 04 


Corrections for Transactions.—In Eq. 37, k vaires from 1 to N not from 
ltoN +1. In Eqs. 46 and 48, the summation is from k = 1 to N, rather than 


fromk =1toN +1. 
333-13 


= 
L; 
6 (91h) 
and 


. 
| _ - 
bss 
ox 4 
> 


DISCUSSION OF STIFFNESS CHARTS FOR GUSSETED 
MEMBERS UNDER AXIAL LOAD 
PROCEEDINGS-SEPARATE 179 


Leonarp K. Stevens.*—The stiffness charts presented by the author 

% Research Asst., Civ. Eng. Dept., Univ. of Melbourne, Melbourne, Victoria, Australia. 
are very useful in the elastic analysis of the behavior of rigid-jointed trusses. 

Mr. Lundquist’s method of stability analysis assumes that the members 
of a rigid-jointed truss remain straight and are centrally loaded up to the point 
of failure, at which point they suddenly buckle. In actuality, the end moments 
introduced by the joint translations resulting from the axial deformations of 
the members produce curvatures, and the strut is not straight at failure. How- 
ever, the ultimate load capacity of an eccentrically loaded, elastic strut is not 
significantly less than when loaded centrally, although large lateral deflections 
may occur before the maximum load is reached. For this reason, it is con- 
sidered that Mr. Lundquist’s method of analysis gives reasonably accurate 
results. Experimental work has shown good agreement in the elastic range. 

It is possible that the secondary stresses produced by the end moments 
may lead to combined stresses approaching the yield point of the material, 
even though the average axial stress is still within the elastic range. If the 
yield point is reached in a slender strut, the load capacity decreases rapidly. 
The determination of the secondary stresses in a stability analysis of a truss 
with struts having medium slenderness ratios is therefore very important— 
especially when the analysis indicates that the average axial stress for limiting 
stability is 70% or more of the yield-point stress. 

Most structures have members in the medium-to-low slenderness-ratio 
range. The yield point of the material is the major factor in determining the 
ultimate loads. As Mr. Goldberg states, plastic flow can sometimes be relied 
on to redistribute the critical stresses. The action when yielding has com- 
menced in a rigid-jointed truss has not been fully investigated. For astrut hav- 
ing a low slenderness ratio, the appearance of initial yielding does not necessarily 
indicate imminent collapse; in fact, the ultimate useful load capacity can be 
considerably higher. An exact secondary-stress analysis can determine the 
load at which initial yielding begins. This load will be the lower limit of the 
ultimate load capacity. The upper limit is fixed by the load producing an 
average stress equal to the yield stress. These limits enable the collapse load 
to be closely bracketed. It is’ in this ‘respect that the charts presented offer 
a valuable aid to the analysis. 
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It is interesting to note the effect of rigid-gusset plates and axial load on the 
moments found in a secondary-stress analysis. Rigid-jointed trusses with 
gusset plates having a ratio of a/L equal to 0.125 and strut slenderness ratios 
of from 55 to 110 have been analyzed by the writer. Near the load at which 
yielding began, the joint moments were as much as 53% greater than the 
moments computed by the classical methods which ignored gusset lengths 
and axial load. 

Joun E. Gotpperc * M. ASCE.—Mr. Stevens’ manifest interest in the 

2} Associate Prof. of Structural Eng., Purdue Univ., Lafayette, Ind. ; 
stiffness charts which have been presented py the writer, and his informative 
and constructive comments are appreciated. His statements in support of 
the contention regarding the effect of gusset lengths and axial loads on the 
magnitude of the secondary stresses are gratifying. It is interesting and 
significant that Mr. Stevens, by using the stiffness charts, has found that the 
effect of gusset length and axial load can increase the secondary bending 
moments by as much as 53%, even in the case of trusses composed of members 
in the intermediate-slenderness range. 

Mr. Stevens is entirely correct in his statement that 


“For a strut (in a rigid-jointed truss) having a low slenderness ratio, 
the appearance of initial yielding does not necessarily indicate imminent 
collapse, * * *. An exact secondary-stress analysis can determine the 
load at which initial yielding begins.” 


A rigidly connected truss is a highly redundant structure, and it appears to 


the writer that one can apply the limit-design concepts and can also invoke the 
general rule that a strut must develop three complete plastic zones before it 
will collapse. 

Corrections for Transactions.—In the first equation of Appendix I, “P,” 
should be “P y’”’ and in Eq. 7, “P”’ (the term preceding the equal sign) should 
be 


~ 
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DISCUSSION OF THE APPLICATION OF HEAVISIDE’S 
STEP-FUNCTION TO BEAM PROBLEMS. 
PROCEEDINGS-SEPARATE 202 


Telemaco van Langendonck.'—It would be more general, to define the step- 
function as: 
Oforx<a 


sen t(x-a) 


1 


lforx >a 
being of no importance the value for x = a, as far as it stays finite (with the 
given formula its value would be 0,5). 

The integral (4) is not general because it is exact only when F (0) = 0. The 
general expression is (c is any fixed value of x): 


Hg f(x) dx = Hg [F(x) - Fla)] + 


The derivative dHa/dx is not equal to zero, as appears in formula (5), un- 
less the considered interval does not include the value x = a. In the same 
way a Hg F(x) only will be equal to Haf(x) if F (a) = 0. 

With these considerations, the application of the step-function can be ex- 
tended to the cases in which there are differential equations to be integrated: 
beams with longitudinal loads, beams on elastic foundations, etc. 

Take, for instance, the case of the beam with two supports of Fig. A, with 
longitudinal load F and transversal load P. The equation of the elastic line is: 


d F M 


where El is the flexural rigidity of the cross section and Mo is the bending 
moment that would exist if there was not the force F: 
Pb 


Mo = Ha P (x-a). (II) 


Putting y' = y" = k?, we have: 


y" + k*y = f(x). 
In general, when f(x) is a polynomial of grade m: 
f(x) = ag x™ + xM-1 4 + + a,x + ao 
the solution of (III) is y = 0, sen kx + C, cos kx + Y, with: 
Y = + xM-1 4... + bx? + + by 


from which 
Y" = m(m-1) byx™-2 + ... + 6b,x + 2b,. 


1. Professor at the Polytechnic School, University of Sdo Paulo, Brazil. 
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Substituting in (II) and equalizing the coefficients of the terms with equal 
exponents, are obtained the coefficients which permit to find Y and, there- 
fore, y: 


a am- 
bm = Dm-1 = bm-2 = 


am-2 - m(m-1)bm (IV) 
eee 


If f (x) contains step-functions 
f(x) = fg (x) + Hg fg(x-a) + Hp fp (x-b) +... 


the integration is made by considering separately each one of the terms 
fo,Hafa, Hafp, ... and by adding the found values of the Y, in order to obtain 
the final Y. 

To integrate 


y" + k*y = Hg[am(x-a)™ + apy (x-a)M-1 4... + a,(x-a) + ag] (V) 


it cannot be used, without modification, the previous process which can be 
applied when there is no Hg in the second member. This happens because the 
derivative of HgF (x) is only equal to zero if F(a) = 0, and this will not happen 
with 


Y = bm (x-a)™ + bm-1 (x-m)™M-1 + ... b, (x-a) + do 


if bo is not zero and, with Y' if b, is not zero. 
To overcome the difficulty it is possible, therefore, substitute b, for 


b,[1-cos k(x-a)] and b,(x-a) for b,[(x-a) - : sen k (x-a)]. So we have 


Y = by)(x-a)™ + ... + b,(x-2)? + b, [(x-a) sen k (x-a)] + by [1- cos 


Y'=mby (x-a)M-1 4 ... + 2b,(x-a) + b, [1 - cos k(x-a)] 


expressions that are equal to zero when x = a, thus having 


d d ' d? d ' uw 
ax (HaY) = Ha & HaY and az (HaY) = Ha Ha Y 


and substituting in (V): 


Hg Y" + Hak?Y¥ = Ha[am (x-a)™ + ... + a, (x-a) + a] 


k?by,(x-a)™ + (x-a)™-! [m(m-1)b,, + +... + 


+ (6b, + k*b,) (x-a) + (2b, + k*b,) = ap,(x-a)™+ ... + a,(x-a) + a 
from which are obtained the same values of bm, bm-}, --. found in (IV). The 
solution wanted for (V) is: 
y = C, sen kx + C, cos kx + H3Y 


with Y given by (VI). 
In the examples of the equations (I) and (II), we have: 


Pb P 


+ k’y = - ET * + H, FI (x-a). 
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and the corresponding to the term Hg za (x-a) is (a, = Hr @_ = 0, from which 


P P 
Di = EM? = Fr Po = 


Ya = [(x-a) -2 sen k (x-a)] 


from which the solution wanted 


Pb P 1 
y = C, senkx + C, cos kx - ** + Ha F [(x-a) - z senk (x-a)} 
where the constants C, and C, are determined with the condition of being 
y = 0 when x = 0 and x= i, that is 


P sen kb 
Ci = Ficsenky 


_ sen kb 


bx H 
- T «= 
y F lk sen k sen kx + Hy (x a) k sen k(x a)| 


which is the equation of the elastic line of the beam of fig. I. The bending 
moments are 


M = Mo + Fy = t (SS sen kx - Hg sen k(x-a)]. 


4 
Pb Pb 
The value of Y corresponding to the term - EI * (a, =- EI’ 20 = 0, from 
which b, = - guy, by = 0): 
Pb Pb 
Yo= - **- 
or 
F F 
Fig.A 
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CORRECTION TO APPLICATIONS OF THE RELAXATION 
TECHNIQUE IN FLUID MECHANICS. 
PROCEEDINGS-SEPARATE 223 


Pages 223-20 and 223-21 should have their numbers interchanged and 
their order of appearance reversed. 


CORRECTIONS TO MECHANICS OF MANIFOLD FLOW 
PROCEEDINGS-SEPARATE 258 


The following is to be substituted for page 258-21: 


The various results are presented as functions of the ratio of lateral dis- 
charge to total and of the ratio of the sizes of the lateral and the conduit. 
Although the results are also dependent on the shape of the conduit and the 
geometry of the lateral, those for circular conduits are thought to be repre- 
sentative of other manifold forms. During a period of nearly half a century 
various experimenters have published data for a variety of other forms, so 
that a considerable amount of information is at hand. In this paper the writer 
has presented new data and systematized the presentation and interpretation 
of the results, It is to be emphasized that tests on a single junction are suf- 
ficient in defining the characteristics of a complete manifold system. For 
many installations, the results presented herein should prove to be directly 


applicable. 
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